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I. INTRODUCTION 


The general aim of this investigation was to study interface stability 
and solute segregation phenomena during crystallization of a model system. 
Emphasis was to be placed on direct observational studies partly because 
this offered the possibility at a later stage of performing related 
experiments under substantially convection-free conditions in the space 
shuttle* 

The study originally eiAvlsaged was a 4- or 5-year program to obtain a 
wide range of careful experimental data in order to test the theories of 
Interface stability. This would have necessitated measurements of a large 
number of material parameters on the system chosen for the investigation, 
the most difficult being the solute diffusion coefficient in the liquid, 
the thermal conductivity of the liquid phase and the solid-liquid 
Interfacial energy. 

In practice, the program was terminated after 3 years which did not 
allow time to complete this detailed Investigation. It has however, been 
possible within the 3-year program to make a number of major advances which 
are important in their own right. The major achievements described in this 
report are: 

1 . The development of a new model system for fundamental studies of 
crystal growth from the melt and the measurement of a range of material 
parameters necessary for comparison of experiment with theory. 

2. The introduction of a new method of measuring segregation coefficient 
using absorption of a laser beam by the liquid phase. 

The comparison of segregation in crystals grown by gradient freezing 
and by pulling from the melt. 
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4. The introduction into the theory of solute segregation of an interface 
field term and comparison with experiment. 

5. The introduction of the interface field term into the theories of 
constitutional supercooling and morphological stability and assessment 
of its Importance. 

The theoretical part of this investigation is presented first. 
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II. THEORY 

A* Interface Field Effects o n Solute Redistribution 

Generally, the solute distribution in a crystal has been described in 
terms of two important parameters, ko which Is called the equilibrium 
distribution coefficient or phase diagram distribution coefficient, and k 
which is called the effective distribution coefficient; ko is the ratio 
of the concentration of solute in the soxid, C$, that that in the liquid^ 
Cl, when equilibrium exists between the two phases at a given 
temperature, while k is the value of this ratio under the actual 
conditions of crystallization. 

From surface studies, one finds that the chemical potential of a mole- 
cule can be quite different in an Interface region compared to the bulk 
phase far from the interface. This is due to differences in local molec- 
ular configurations between the two domains arising from (i) gross struc- 
tural changes in the molecular potential function for an atom situated 
there rather than in either bulk phase, and (ii) inhomogeneous fields of 
stress, electrostatic potential, magnetic potential, etc. Therefore, in an 
interface region, the Gibbs free energy is modified by introducing an 
"Interface Potential" term through an extended chemical potential 

nJ _ pj 

= + kT £n a^(x) + AG°^(x) , (1b) 

where j refers to the solute, a is the activity, p® is the standard 

state chemical potential, AG°^(x) is the interface potential which is the 

t/hi 

spatial variation of the energy of interaction between the j solute and 
all of the interface fields and k is the Boltzmann constant. 


« • 
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In an Interface region, the chemical potential difference which comes 
from differences In local molecular configurations Is called the "Interface 
Potential: and this produces the "Interface Field," The major Interface 
fields of Interest are categorized by writing 

SG°(x) = f(g(x), o(x), H(x), <|)(x)) . (2) 

Here g(x) accounts for the variation of local order, local coordination 
and local density with distance: cf(x) accounts for the variation of local 
stress; H(x) accounts for the variation of local magnetic potential; and 
<j)(x) accounts for the variation of local electrical potential. At this 
moment, the formulae of all Interface potential terms are not well 
defined. However, the interface field contribution can be written 
explicitly [1] as 

AG°(x) [| (V^)2 + 5^] + + ... ( 3 ) 

where Z Is the valence, <f» is the electrostatic potential, e is the 
permittivity, a is the stress, E* is the elastic modulus, and Aph.M. 

Is the chemical potential change relating to different bonding energies for 
solute versus solvent species in the Interface region, and neglected 
contributions include activity coefficient effects. 

In order to evaluate the effect of an interface field on the solute 
distribution, consider a frame of reference moving at constant velocity V 
with the Interface subject to a nonstlrred or a well-stirred fluid at 
distances greater than x = 6 from the interface and subject to interface 
interaction energy contributions of the form 
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6Gg(x) s 0g exp(cgx) , 

X "< 0 , 

(4a) 

6G°(x) s exp(-o^x) , 

X > 0 . 

(4b) 


This fom of field decay has been chosen for Its simplicity, where p and 
<9 are constants. The overall solute distribution might be represented 
qualitatively as In Fig. 1 for a stirred fluid. The flux of solute, 3 is 
given by 

3 S -d3c - ~ ^6G°(x) - , (5) 

kT 

1{> a coordinate system moving with the Interface at velocity V, In the 
steady state, 

3 • 3 = 0 , (6) 

which has the general solution 

VC + + Be“'^)c = X , (7) 

where 

for solid: B s - o,' = a (8a) 

kT S 


for liquid; 


(8b) 



Figure 1 Schematic representation of solute distribution in the solid 
and^liquid during freezing from a convectively stirred liqui 
fx > 6) and with an interface field for the case of k < 1 
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Using the integrating factor technique, Eq, (7) has solutions of the 

form 

C s E" cxp(- ™ x) exp(- "It + A exp(- | x) cxpC- "It 

X / exp(~ t) exp(~, e“'*^)dt . (9) 

Solutions for C$ and C|_ are plotted in Fig. 2 as a function of 
interface position x for a range of growth velocities and with parameters 
tabulated in the figure. For details of the calculation, see references 
t1] and [23. 

The logarithmic scale is necessary to reveal the two very different 
scales of phenomena operating in the liquid; l.e. , the phase diagram solute 
partitioning profile and the field-affected interface segregation profile. 
As a first approximation, we might think of the two as acting Independently 
and Fig. 2 is a direct superposition of the two effects. Thus, the solute- 
depleted zone at small x for large V indicates a measure of the time 
response of the solute ahead of the Interface to the interface field. At 
higher velocities, a shorter time exists for such a response so the mlnumum 
moves closer to the interface. Only when V ” 10"2 cm/sec does any 
appreciable pile~up occur at the interface and by V “ 10“^ cm/sec, the 
equilibrium pile-up is developed. The solute distribution in the solid 

_i 

behaves as expected with the profile extending a few Og and with the 
solute trapping effect becoming complete at V “10"”*. 
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We seek an expression for the effective distribution coefficient, k, glv<'n 
by 


k P Cg.-AXg)/CL(«) . 


(10) 


for a stirred melt such that Cj^ p x 6. The solution [1,2] is 

of the form 


and an approxifnate summation for the gamma function Y5 gives 

k. 


(11a) 


k = 


k + [expf— 
i kT 


1 <1*V - k K-^ ) 

S 1 tlOu 6 l 


iJt ,-V/(0LDt) 


\1 i « 

(11b) 


In the limit as $5 and P(_ +0 , ys ’♦’ 0, 

V kT 


f. ill 1 


Yl -^1 ’ 


and ki ■*■ ko» Thus, we have 


k = 


ko + (1-k^) e 


-V6/D, 


(12) 
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whic5i Is the familiar BPS result [33 for no Interface field effects. In 
Eq. (11), Is the Interface distribution coefficient Cg(0)/Cj^(0) which 
can be obtained In terms of by putting Pg = nj_ (see Eq. (1)). An 
alternative form of Eq. (11b) is 


k = 




k' f- -X_ 


)(- 


kT 

P. 


) YlJ exp{- 


V5, 


(13) 


where kj^ = kQ[Yg(“)/Y5(0)3 [Y|_(“)/Y|_(0)]"^ . 

To illustrate the effect of 3, a and D changes on k, we first 

explored the range -5 < (3/kT) < 5, 10^ < < 10^, 10”^ < V < '0"^ with 

Dg = 10"®, Dl = 10"®, k^!) := 10‘\ 6 = 10'”' and 10^ < < 10^. If 

k'/cs, D, « l/acD-, and 3- effects are dominant over the 3, effects and 

vice versa If » 1/ogDg. As expected, for 3g < 0. ;;.v3/kT) > k(0) 

while for 3g > 0, k(3/kT) < k(0). This effect is illustrated in Fig. 3, 

which plots k(5), k(0) and k(-5) as a function of V for 3|_ = 0, 

6 2 

Og = 10 and = 10 . Only a small difference occurs between k(-5) 

y 

and k(0). Usually, k + k^ as V 10" cm/ sec for all 3$ and 05 
k 1 as V 10"® - 10"^ cm/sec. 

There does not appear to be any suitable data with which to 
quantitatively compare our theoretical results primarily because so much 
parameter Information is needed to meaningfully test Eq. (13). Only with 
this objective in mind would investigators gather and record all the needed 
data for an adequate test. However, it is useful to pro'^ect how one would 
go through such a comparison with experimental data and we shall do this 
with Hall's data on germanium [4]. To test Eq. (13), we would plot 
Zn (1/k - 1) versus 6 for fixed V and versus V for fixed 6. In the 
former case this should be a linear plot with slope V/Dl and intercept 
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(l+Yg) 

^n[— rA 


W ir “V/(a,D, ) 

(•„-Vk- f \i • 


'W' 


('!*) 


This Interesept should vary with V and the magnitude of the variation will 
indicate the magnitude of the field effect. In the latter case a nonlinear 
plot should obtain for a nonzero field effect and the greater the degree of 
nonlinearity^ the greater the magnitude of the field effect. 

The simplifying assumption for Ge or Si systems is that, because 
the liquid is metallic, we can approximate ~ 0 so that 


1+Yc 


«n(l- 1) . 1) - ^ . 


(15) 


In addition, assuming that is approximately independent of V, 


dV 


+ (_1 ) (_I§) 

), ^ M-k'^ ^dV N +0 


. 1 « Pc/kT ^ 

D| 1-k' ^oicDcn n! J 

L 0 n=i S S 


(16a) 


6 




, for |Ps/kT| « 1 (16b) 


so that the initial slope yields a value for Pg/og through its departure 
from the zero interface field value of -6 /Dl. Finally, as V 


dYg/dV 0 and 


[d An(-j^ - 1)] 
dV 


6_ 

D. 


(16c) 


Ki cmmC. 
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Replottlng Hall's data In this fashion, as In Fig. we find that 
d tn(1/k - 1)/dV seems to approach a constant slope at large V. From 
this value, we obtain 6 /Di_ " 57 sec/cm and this yields a value of 
~ 2.4 X 10"^ cm^/sec for Sb in Ge. From the V - 1 inch/hr 
intercept, we obtain " 2.2 x 10" compared to a tabulated value [5] of 
k^ = 3.0 X 10 " . From the slope at V = 0 (same value as V = 1 In/hr), 
which is -428 sec/cm, we obtain ( gg/kT) (1/o^Dg) w -1.7 x 10^ sec/cm by 
assuming that ( 35 /kTl « 1. From Fig. 4 at V = 'J in/hr, we require that 
Jln( d+Yc-k' )/k' ) be "5 so that Yc " -0.67. ‘fhis leads to 
Og = 9.2 X 10^ cm'^ (Dg = 5 X 10"^^ cm^/sec) [ 6 ] and p 5 /kT •> -1.17 for 
the ( 100 ) face, which violates our initial assumption of |pg/kT| « 1. A 
more correct analysis yields values of Sg/kT -1.7 and ®^^g •> 2.62 x 
10 “^ cm/sec for the ( 100 ). 

The values of Bg/kT for the various crystal faces differ by "0.1 
with the (111) having the largest magnitude. The sign of the binding 
energy indicates an Sb buildup at the Interface probably due to a 
negative electrostatic potential binding of Sb+ ion cores to the 
interface region. 

If Pl is not neglected, it should be possible to separate the 
effects of interface fields in the solid and liquid phases but it was 
difficult to find the most suitable values for gg, g^, CgDg and 


using the data and computer facilities available. 


ti (k 



Figure 4: Comparison of 



15 - 


B , Interface Field Effects on Constitutional Supercooling 

Generally, the constitutional supercooling criterion (CSC) may be 
written as 





x=0 


( 17 ) 


where Is the liquidus temperature and is the actual temperature. 

For the steady state growth situation with no Interface field present, this 
CSC becomes 





°L '<0 


(18) 


Thus far, this CSC has been based on a static analysis, provided that 
Interface attachment kinetics are infinitely fast and the Interfacial free 
energy is negligible. However, when we have an Interface field present, 
Eq. (18) is no longer valid and we must modify the liquidus slope to 
which includes the interface field effect: 


M* = exp(6G°/RT) , (19) 

where 6G° = (AG^^ - CAG^^), and superscripts 1 and 2 refer to solvent and 
solute, C = v°^/v°^ is the volume ratio for the liquid phase and v° is 
the specific volume. 

Using the same interface field representation defined above, and noting 
that C. = C-/k. , we find [2,7] that 

X ^ X 


(dT/dx). 


-M* 


V - D. 


exp(- 

Cj— [1 


[1 

kT 


00 

- I 

n=1 


(Pg/kT)'^ 


(1 




■] - 1 } 


)n! 


( 20 ) 


For the case of metallic liquids, »• 0 and o Is large so we 
probably need to keep the term in Eq. (20) even though we may be 

justified in setting exp(-B|_/kT) " 1. For Bg = 0, this would tend to 
reduce the r.h.s. of Eq, (20) and thus decrease the tendency toward 
constitutional supercooling. Depending on the sign and magnitude of Bg, 
the tendency can be either further decreased or increased. For the case of 
germanium discussed earlier, since Bg/kT " -1.3 and ko' ** 0.73 ko, w 
expect that, for B|_ = 0, the r.h.s. of Eq. (20) will be decreased below 
that of Eq. (18) so that the Interface should be stable even in the 
presence of some constitutional supercooling as evaluated using the old 
criterion. Estimates of the size of this effect for CsCdCla under various 
assumed values for Bg and Bj_ (see Section VI) are shown in Fig. 5. 


C. Interface field Effects in the Morphological Stability Analysis 
An alternative and elegant approach to the planar interface insta- 
bility criterion was developed by Mullins and Sekerka in 1964 and indepen- 
dency by Voronkov [9]. This treats the time-dependence of an infinitesimal 
distortion of a plane interface. If it grows, it causes a permanent shape 
change and leads to an unstable interface state; whereas, if it shrinks, it 
leads to a stable interface state. Mullins and Sekerka assumed a 
sinusoidal disturbance of amplitude 6 and calculated the conditions under 
which 6 (= d6/dt) was positive. 

We consider the solute redistribution for both solid and liquid in a 
frame of reference moving at constant velocity V\ with the interface 
subject to a nonstirred liquid and subject to an interface Interaction 
energy contribution of the form 





- * * 
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5Gg(Y) » Pg exp(cgr) , 

r < 0 

(21a) 

6G°(y) s 3j_ expNa^^r) , 

r > 0 

(21b) 


where 3 and a are constants and subject to a sinusoidal shape 
perturbation with an amplitude 6 and wavelength X. 

The governing equation for solute transport Is 

V 

5 = -D[^C + T- + — C] (22a) 

solute D 


where D Is the diffusion coefficient of solute, which Is assumed 
constant. The governing equation for heat transport Is 


5 


thermal 



(22b) 


where Dt Is thermal dlffuslvity which is also assumed to be constant. 

A steady state condition will be assumed and this requires that, 

ac/8t = -^‘^^soiute ' ° ^ "^‘^^thermal "" definition 

of 6G° and only taking the first order 6 terms into consideration 

?(SG°) = (50°) 

= aP - aS cos wZ) + i<(w5 sin iiiZ)) , (23a) 

and 

^•^(A°G) = o3 e“^{a + o)^ - oi^)6 cos uZ} . 




(23b) 
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Unfortunately, using these formulae, we cannot solve analytically the 
governing equation of solute transport which has the ultimate form of 


,2c t 

kT 


+ ;r ‘ 

kT D 3X 


a 0 


( 24 ) 


However, once we specify a certain svitem, we can solve the above equation 
by numerical analysis methods. 

To solve the solute transport equation analytically, we assume that the 
Interface Interaction energy has a linear form such as 


6G° = 0(1 + ar^) , (25) 

and this holds only within |r^j < |1/“|* Then, we can write 

^(6G°) s 0o(t + i< 6u) sin wZ) , (26a) 

^•^(6G°) = 6a6a)^ cos ml , (26b) 

and we have the general solution to Eq, (24), which is 

= Aq + exp[-(— + -~)X] + (A 2 + A^ e’“**^) cos ml . (27) 

kT 

If we obtain a similar solution for |x| > jl/aj and apply continuity 
equations for C and 3 at x = jl/aj, we can obtain concentration and 
temperature distribution; After much tedious algebra [2,7], the stability 
relation becomes 


41 « 
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4>a-t “ <Sf> * 1<?> 

''x 

C„(^) (NUM) 

+ to(Kg+K|^) ” M* — •} ] * ( 28 ) 

where (NUM) and (DEN) are complex functions of U5» <*3* V etc. To 
Illustrate the Interface field effect on 5/6. we will identify the Mullins 
and Sekerka result as (6/6)j^_g and ours as (fi/S)^_j* We first explore 
the range -3 < 3g/kT < 3 and -3 < 3^^/kT < -1 x 10’^ for 10"^ cm/sec < 

< 1 cm/sec with cxg = Oj^ = 10^, and w « 10^ for the Ice-water-KCl 
system. We used all the experimental values for the Ice-water-KCl system 
quoted by Kotler [10] except for the surface energy which is y s 29.1 m3/m^ 
from Trlvedl [11]. The values used are AH = 80 cal/g, T^ = 273°K, 

Ks = 4 Kl = 5.3 X 10"^ cai; (cm/sec - “C), k© 5 x 10"**, 

Ml = 1 .7“C/(mole-llter"M » Coo = 0.01 mole/llter, Dl = 5 x 10"® cm^/sec, 

Ds = 10"® cm^/sec, astd Gl s 25'’C/cm. From Fig. 6 we can observe the very 
marked difference between the two criteria. This difference is not easily 
explainable by using a single special instability mechanism but it should be 
obvious that the interface field effects influencing the constitutional 
supercooling criterion, and the point effect on the lateral diffusion, are 

f 

two major contributions of this very large difference between (6/6)/^_j and 
(5/6 )|v|_ 5. From Fig. 6, we can especially note that only a small Pl/^T 
value is needed to generate a large effect. Therefore, dealing with a 
material system which obviously manifests an interface field (electrostatic), 
but neglecting the field effect on the stability criterion as people have in 
(he past, produces a markedly erroneous result. 



Figure 6; Calculated plots of versus V illustratin, 

the importance of (3g/kT and B^/kT for Ice-Water-^KCl system 
and calculated plot of (<5/6)j^_g vei'sus V for this system 
C3g/kT = 0, B^AT = -1 X 10'4) 
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III* CHOICE OF MATERIAL FOR EXPERIMENTAL STUDY 

The material requirements sought were that crystal growth should be 
possible by the Czochralski method and by gradient freezing, ‘ 
transparent crystals could be grown without major problems and thiSt the 
melting point should be sufficiently low tc permit direct observ/jtlon of 
the growing crystals. In addition, it was necessary to find a suitable 
dopant of fairly low phase diagram distribution coefficient for 
quantitative studies of solute partitioning during Czochralski growth. 

These requirements were met by cesium codmium chloride, CsCdCla and a 
preliminary report on the favorable properties of this material for 
crystallzatlon studies has been published [12], 

A. Phase Diagram 

The phase diagram [13] for CsCl - CdCl 2 (Fig. 7) shows three compounds: 
CsCdClg which melts congruently at 553®C, Cs 2 CdCl 3 with a congruent 
melting point of 488"C and CsaCdCls which melts incongruently at 390'’C. It 
has been confirmed by DTA that CsCdCla melts congruently with no phase 
transitions between the melting point and room temperature but the melting 
point of purified material is actually 545“C. 

Investigation of the crystal structure of CsCdCla has confirmed that it 
is hexagonal with a = 7.418 A and c = 18.39 A as reported by Siegel and 
Gebert [14] rather than monoclinic as given by Naray-Szabo [15] and 
Belyaev, et al, [13] or cubic perovskite as reported by Marsh and Savage 
[16], Molten CsCdCla has appreciable volatility and a quartz crucible, 

42 mm in diameter and containing about 160g of material at its melting 
point, was found to lose between 0.3 and 2.6 wt%/day depending on its 



deviations from stoichiometry (when we used the materials as received from 
vendors). The CdCiz-rlch material of highest volatility in this range had 
a melting point of 535*C and lost CdCla by evaporation, but the 
stoichiometric material was found to evaporate congurently. The solubility 
of CsCsCla in water is reported as 53.7 g/Jt at room temperature [16], 
Although this rather high solubility suggests that powdered material will 
be hygroscopic, the rate of attack by atmospheric moisture of the surface 
of crystals was found to be very slow and no signficant handling problems 
were encountered. 

The source materials used were CsCl of 99.95% purity (Harshaw) and 
"ultra pure" CdCl 2 • 2 H 2 O (Alfa Ventron). The cadmium chloride hydrate was 
dehydrated first by a very slow vacuum treatmefit at room temperature to 
prevent hydrolysis \CdCl 2 + 2 H 2 Cd(0H)2 + 2HC1 ) and then by heating in 

vacuum at 200“C for 24 hours. The dehydrated CdCi 2 was purified by zone 
refining or distillation (vapor transport). Both purification methods were 
applied to CsCl but the vacuum distillation gave higher purity material and 
was, therefore, used in a majority of the experiments. Major impurities in 
the distilled material were 1 ppm Ca and 0.4 ppm Si. 

Problems were encountered in zone refining CsCl because of the stress 
generated at the solid-liquid interface during the 2nd pass which 
frequently led to fracture of the silica tube. This problem was alleviated 
but not eliminated by coating the inside of the tube with pyrolytic 
graphite. Vacuum distilled CdCl 2 contained 3 ppm of Ca, 1 ppm of A1 and 
< 1 ppm of Mn and Si. Some experiments were performed using a commercial 
source of high purity anhydrous CsCl and CdCl 2 (Metalsmart) . Impurities 
reported for these materials were 2 ppm Si and 0.5 ppm Cu and Mn in the 
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CdCIa and 4 ppm Ba and Rb, and 1 ppm each of Al, Ca, Cu, K, Mg, Na and Si 
In the CsCl. 

B. Solute (Dopant) 

CsCdCla can be doped for solute redistribution studies by replacing 
CdCl 2 with a divalent transition metal chloride. C 0 CI 2 produces a blue 
coloration, CUCI 2 an orange-brown and NICI 2 a pink. It is assumed that the 
transition metals enter the lattice as divalent ions but this has not been 
confirmed. All these impurities exhibit a low segregation coefficient but 
quantitative data are not available except for Co. The rate of back- 
diffusion of copper ions during cooling of the crystal may be very fast, 
depending on the orientation. Therefore Cu is not a suitable dopant for 
this study. Absorption spectra of doped crystals in the visible and near 
UV regions shows, in each case, strong absorption at short wavelengths 
(approximately 3000 A). In the case of the Cu-doped crystal, a broad 
absorption band extended to about 5000 A, with no other peaks in the 
visible range. The Co-doped and Ni-doped crystals had peaks centered at 
6200 A and 5000 A respectively. Since the absorption peak of Co-doped 
CsCdCla Includes the 6328 A emission line of the He-Ne laser and cobalt 
does not exhibit appreciable back-diffusion, this dopant has been used for 
this study. 

C. Property measurements 

We attempted to measure all the experimental parameters for the CsCdCls 
system needed to compare the theory with the experimental results for 
interface instability of CsCdi.^Co^Cla . Most of the parameters were 


• <i • 
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found successfully? however, insurmountable difficulties were found in 
measuring the thermal dlffusivlty. The parameters that were successfully 
measured are listed below. 

1 . Llquidus slope measurement 

The slope of the llquidus curve was measured on a Dupont Differential 
Scanning Calorimeter (DSC). Great care was necessary in sample preparation 
in order to eliminate moisture and to produce homogeneous samples of known 
composition. Because of the hi^ih degree of supercooling encountered wJth 
this material, the llquidus temperature data were taken during heating of 
the sample from the onset of melting. The results are shown in Fig. 8. 

2. Specific heat Cp 

The specific heat of CsCdCls was also measured by using the Dupont DSC, 
with sapphire as reference. The specific heat was found to be approxi- 
mately 0A5 3g“^ deg"^ around room temperature, 0.33 3g"^ deg“^ just below 
the melting point, and 0.30 ± 0.01 3g"^ deg"^ in the temperature range from 
580 to 620°C. 

3. Latent heat of fusion AHf 

The latent heat of fusion of undoped CsCdCls was measured by the DSC 
and found to be 132 3/g. This experiment was done simultaneously with the 
luquidus slope measurement. The latent heat of fusion of CsCd-i.^Co^^Cls 
showed a slight Increase with Co concentration. 
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4. Density p 

The density of molten CsCdCls was measured by a direct method using a 
cathetometer to measure the height of an Isothermal column of liquid. The 
density was found to be 3.03 g/cm at 564 ± 4“C. The error of measurement 
of the density Is about 2%. 

The density of solid CsCdCla was determined at room temperature by an 
Archimedean method. The value obtained was 4.03 g/cm . It was found that 
the GsCdCla expands 24.7% of its volume when It melts, which is rather 
large In comparison with most molten salts. 

5. Viscosity n 

The viscosity of molten CsCdCla was found to be too low for the falling 
sphere method, so an Inclined tube technique [17] was used. Using water at 
room temperature as a standard, the viscosity measured using bronze spheres 
was found to be n = 1.35 cp just abn"® the melting point. 

1 

6. Thermal diffusivity Dt 5 

1 

We tried to measure the thermal conductivity of liquid and solid j 

■ i 

CsCdCla using a "semi-infinite plate" method with periodic temperature j 

i 

variations. The liquid (or solid) specimen is contained in a double : 

crucible with a relatively thin layer of material between the upper and 

lower crucibles. A periodically varying heat source supplied heat to one s 

surface and the temperature was measured by axial thermocouples welded to 

the two metallic surfaces which enclosed the sample. The apparatus was j 

i 

i 

calibrated against standard material which was sodium nitrate. However, | 


■■r NHar.. 

m 









when we experimented with CsCdCia in Ni crucibles, it severely attacked the 
crucibles and we could not find suitable crucibles for CsCdCla. 

The apparatus was reconstructed using silver, but CsCdCla vapor was 
found to attack the heating element very rapidly once the temperature was 
raised above the melting point. Time was not available to build a modified 
version with a protected element, and att-?mpts to locate a laboratory with 
the capability to undertake this measurement were unsuccessful. 

Table 1 lists the measured parameters for CsCdCla. 
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Table 1 

Measured Parameters of CsCdCla 


Liquidus slope 


6/“C/wt% 


Cp Specific heat 0.45 3/g-“C at R.T. 0.30 ± 0.01 3/g-“C at 600“C 


132 3/g 


4.3 g/cm^ (solid) 3.03 g/cm^ at M.P. (liquid) | 


1.35 cp at M.P. 



AH Heat of fusion 


p Density 


U Viscosity 


IV. STUDIES OF INTERFACE SHAPE AND STABILITY 


Studies of cobalt-segregation and the onset of Interface Instability 
have been made using the apparatus shown in Fig. 8. A 2-zone furnace was 
wound on a grooved ceramic tube which had few vertical slots 7 cm long to 
allow illumination and observation of the growing crystal. The temperature 
gradient was controlled by two independent temperature controllers, and 
crystal growth was achieved by programmed cooling of the upper zone. 

Theories of solute distribution and interface stability normally 
apply to a plane Interface and so attempts were made to achieve a fairly 
planar solid-liquid interface over a wide range of values of g, the 
fraction of metal solidified. The general principles which determine the 
interface shape are well established [18]. If the interface is to be 
planar, the heat flux from the liquid must be the same as heat through the 
solid; if the crucible receives a net heat flux from the exterior, the 
interface will tend to be convex and if it radiates heat to the surround- 
ings the interface will be concave. Solution of the heat-transfer equation 
has been given a sample of infinite length [19] but only a very approximate 
solution is possible for a finite length. In practice it is found 
difficult to maintain a heat interface as growth proceeds [19] and typical 
interface shapes at various stages of growth in unseeded and [0001] seeded 
cases are shown in Fig. 9. This shows typical trends in the development of 
the interface shape of CsCdCls with time as growth proceeds at slow speed 
(< 1 mm/hr) in a temperature gradient of 22"C cm~^ in samples containing a 
low concentration of cobalt (” 0.005 a/o Co replacing Cd). The interface 
is normally slightly convex (as seen from the liquid) and the convexity 
tends to decrease as solifification approaches 100%. In the case of the 
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Figure 9 . General view of furnace and 
controls. 
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seeded growth (Fig# 9b) the convexity passes through a maximum In the 
region where the axial temperature gradient In the liquid Is a minimum as 
suggested by growth rate data. The convex Interface shape suggests that 
Ks/K-| > 1, although this has not yet been confirmed by direct measure- 
ment. 

Our observations have shown that the convexity of the Interface 
decreases as. the furnace temperature gradient Is Increased. In the furnace 
used to date, the maximum temperature gradient is H2.2*C cm"\ and this was 
found to be adequate to give a substantially flat Interface over a wide 
range of growth conditions. Reducing the gradient to approximately 14"C 
cm"^ caused pronounced convexity of the interface under similar condi- 
tions. The dependence of interface shape on growth speed was relatively 
minor, flatter interfaces being observed as the cooling rate was Increased 
to the range where the interface shape became unstable and cellular growth 
occurred. 

The most striking observation in our preliminary experiments was the 
occurrence of a marked dome shape as the cobalt concentration was increased 
from 0.005 to 0.01 a/o (see Fig. 10). At such a low concentration the 
cobalt cannot be expected to cause a sttong change In any of the crystal or 
melt properties. However, the cobalt does cause a very noticeable Increase 
in the coloration of the melt and therefore in the absorption coefficient 
over a range of visible and presumably Infrared wavelengths, and this 
increased absorption increases as the crystal grows and cobalt is rejected 
into the liquid by the normal segregation mechanism. The occurrence of a 
markedly convex interface as the heat flow through the liquid phase is 
decreased is entirely consistent with theory: a diminished axial flow is 
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Figure 10 (a) s Variation in interface shape of CsCdCl^ as 

growth proceeds at slow speed in a temperature 
gradient of 22.2°C/cra in samples of low cobalt 
concentration 005 at % Co replacing Cd) 

(a) unseeded growth '(.b) JOOOl] seeded growth. 
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compensated by increased radial flow of heat into the sample under 
conditions where the flow through the solid remains constant, so that the 
convexity increases. 

It was reasoned that decreasing the axial heat flow through the 
crystal would oppose this effect and therefore restore a more nearly planar 
Interface. An experiment was therefore performed in which the bottom 
surface of the capsule was coated with a gold film to reflect Infrared (and 
visible) radiation transmitted through the crystal., In contrast to our 
predictions, the presence of this gold coating caused an even steeper 
convex interface, and the crystal grew with a liquid boundary groove 
separating it from the capsule along its whole length (see Fig. 11). The 
reason for this effect is not understood, although conduction of heat 
through the silica crucible may be an Important effect in maintaining a 
relatively high temperature at the walls. 

The growth rate during growth at constant cooling rate was found to 
vary slgnflcantly during the growth of a crystal typically 5 ± 1 cm in 
length. Typical data is shown in Fig. 12. In the steady state, the growth 
rate V is given by 


V = 0/G (29) 

where 0 is the cooling rate and G the temperature gradient. Although 
this equation is not expected to apply to a relatively short crystal, a 
maximum in the growth rate is to be expected since G will Increase at 
both ends of the sample due to radiation from the hot zone and to the cold 
zone of the furnace respectively. In addition, the steeper decrease in V 
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Figure 11: 


Change in interface shape of CsCd« onCo^ «tC1- 

u • y 7 U)» ui, j 

at same stage of growth (a) original condi- 
tions, see right hand figure of Figure 10(b); 
(b) with a gold coating across bottom of 
crucible. In this latter case the crystal 
grew surrounded by a groove containing liquid. 



DISTANCE ALONG CRYSTAL (mm) 


Figure 12: Variation of grov th rate of CsCdClo crystals 

with position during growth by gradient 
freezing. Temperature gradient: 17.9°C/cm; 

0 22.2°G/cm. 







during the final stages of solification Is due to the bulld-up of cobalt 
concentration ahead of the interface, with corresponding decrease in 
liquidus temperature. 

No attempt was made to explain quantitatively the variation of the 
growth rate as solification proceeded. Both this and especially the 
changes in interface shape are considered worthy of more detailed investi- 
gation, which would involve measurements of the temperature profile within 
the melt. The Influence of heat absorption in the liquid phase as the 
cobalt concentration increased is of particular interest. 

Table 2 summarizes our initial investigation which was aimed at 
identifying the critical growth rate for the onset of instability and its 
dependence on the temperature gradient to which the sample was subjected. 
The same [0001] seeded sample was used in each case, with the cobalt con- 
centration maintained constant at x = .0025. The results showed a general 
tendency for the maximum stable growth rate to increase with the tempera- 
ture gradient, as expected. The "approximate growth rate" values quoted in 
the table are unreliable as they Include the region where growth slows 
signficantly towards the end of the solification. If we take the critical 
cooling rate in each case and assume Eq. (29) for the growth rate, the 
maximum stable growth rate is 2.2 mm/hr for a gradient of 22.2°C/cm, 
mm/hr for a gradient of 17.9“C/cm and 1.0 mm/hr for 14.3°C/cm. Hence 
Vmax “ approximately, in agreement with theory. 

These early studies showed, however, the need for continuous monitor- 
ing of the interface position in order to obaln an accurate value for the 
growth rate at the instant of breakdown. It was found to be difficult to 
determine precisely the moment where growth became unstable, even on film. 


Table 2 

Summary of Obsrvatlons of [0001] Crystal 


A. Temperature Gradient 22.2°C/cm 


Cooling Rate 
(“C/hr) 

Appro>:imate 
Growth Rate 
mm/hr 

Range of 
Measurement (mm) 

Observations 

0.8 

0.65 

8-22 

Stable 

2.5 

0.9 

31-34 

stable 


1.6 

14-37 

Perhaps temporary 
instability 

4.5 

1.3 

19-32 

In verge of 
instability 

4.75 

1.4 

23-34 

On verge of 
instability 

5.0 

1.8 

19-35 

Slight Instability 

Temperature Gradient 17.9‘’C/cm 



1.0 

0.5 

9.30 

Stable 

2.0 

1.0 

15-35 

Stable 

2.25 

1.0 

24-35 

Stable 

2.5 

1 .2 

23-32 

Slight instability 

4.5 

1.8 

26-32 

Unstable 

Temperature Gr 

adient 14.3°C/cm 



0.66 

0.6 

16-24 

Stable 

1.0 

0.7 

24-31 

Stable 

1.25 

0.7 

20-32 

Stable 

1.5 

1.0 

20-27 

Unstable 

2.0 

1.0 

15-18 

Unstable 
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More careful observations suggested that the above values for 
low by a factor of 2. In addition > the validity of the results at the 
lower gradient was affected by the curvature of the Interface which became 
pronounced when values below 22" C/cm were used. 

Since the constitutional supercooling criterion is of the form 



mk(1-kp)C 

■ 


m(1-k^)C„ 


(18) 


a test of l^'s validity requires measurement not only of G and V^ax 
but also of Coo (in addition to m, kg and D|_ in separate 
experiments). In all previous experiments of which we are aware, Coo has 
never been measured directly during the course of a crystal growth 
experiment but has been inferred from later measurements of the dopant 
concentration Cs in the crystal. Emphasis was therefore given to the 
development of a method for the determination of the concentration Cl 
of cobalt in liquid CsCdCls since this concentration can then be measured 
on a continuous basis. The same method can be used for measurement of the 
distribution coefficient k, Including (by extrapolation) the equilibrium 
value ko- 

Refinement of our filming technique continued in order to improve the 
visibility of the onset of Interface breakdown, and seed crystals which 
would completey fill m 8 mm I.D., fused silica tube of square cross 
section were Introduced. The use of high magnification confirmed that, 
when the interface does become unstable, it breaks down into the cellular 
structure predicted by morphological stability theory [8], An example is 
shown in Fig. 13. 
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Figure 13: 


Solid-liquid Interface during crystal growth of 
CsCdCl- showing onset of cellular growth. The 
periodicity of the cellular structure is about 
200 um. 


V. MEASUREMENT OF DISTRIBUTION COEFFICIENT USING LASER LIGHT ABSORPTION 


It was mentioned in Section III.B that cobalt-doping of CsCdCla intro- 
duces a broad absorption peak centered at 6200A and that this peak includes 
the 6328A emission line of a He-Ne laser. This suggested that the absorp- 
tion of a He-^Ne laser beam could be used as the basis of Co determination 
in molten (or solid) CsCdCla and an arrangement was set up using a Hughes 
10 mW laser (Model 322SH-PC) and a detector (Model SD-W4111261) from 
Silicon Detector Corp. The output from the detector used in a photovoltaic 
mode was found to depend very strongly on the cobalt concentration. 

The absorption was calibrated using molten Cs Cdi^x Co^ Cl 3 
samples of known values of x, and the absorption coefficient a was found 
to depend linearly on the Co concentration. The voltage registered by the 
dectector was around 5V for x = 0 and fell to ~ 1 mV for x = 0.004. 

The measurement is therefore very sensitive within this region but the 
signal- to-nolse ratio fell below the useful range for x > .004. 

A minor disadvantage of the laser absorption method is that the 
absorption coefficient for the He-Ne laser light is temperature-dependent. 
Since this dependence was found to be significant, the calibration was 
carried out at a temperature just above the melting point. A second 
calibration was obtained at ~ 630“ C, v^hich decreases the apparent value of 
a by 10%. This decrease is attributed to broadening of the absorption 
peak due to motion in the fluid. 

In order to obtain accurate values of distribution coefficient, it was 
found to be necessary to measure the length of both the liquid and solid 
column using a cathetometsr. The ratio of these two quantities gives g, 
the volume fraction of material solidified. 
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The absorption coefficient In the liquid just ahead of the interface 
was used to give a colbalt concentration Cl as a function of g. 

According to the normal freezing equation, 

Cl = S^d-g)'^"'’ , (30) 

where Sq Is the amount of solute In the liquid at the onset of growth 
and k Is the effective segregation coefficient. Our data was therefore 
plotted as graphs of In Cl versus In(l-g) which have a slope of 
-jk-l|. A computer least squares fit was used to find the value of k for 
different values of the growth rate V. Sample data were shown In the 
report for November 1981. 

Data over a side range of values of the growth rate were obtained for 

both <0001> and <34?3> orientations. The latter was not chosen 
deliberately but measurements were made on a crystal which grew after a 
<0001> seed had melted due to a sudden Increase In temperature. 

Figure 14 is a summary of the experimental data In the form of plots \ 

of An(k"^-1) versus growth rate V for the two orientations. This form ■ 

of extrapolation Is suggested by the BPS equation [3] 



where ko Is expected to be the equilibrium value of the distribution 
coefficient. Although Eq. (12) Is not expected to apply exactly to 
gradient freeze growth but fits the data within the experimental error. 
The values obtained are ko = .02(7) for <0001> and ko = .07(4) for 


<3473>. The estimated error is ± .02 for <0001> and ± .03 for <34?3> so 
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the difference between the two directions Is Just outside the estimated 
limits. 

The slope of the least squares fit Is 2.0 x 10^ s cm"* for <0001> 

and 2.1 X 10^ s cm"* for <34?3>. This agreement is fortuitous but the 
use of D|_ = 7.9 x 10"® cm^ s"* (see Section VI) gives 6 = 0.16 mm which 
is probably correct to the order of magnitude. 

Attempts to detect a solute boundary layer of high Co concentration 
ahead of the Interface were unsuccessful as would be expected if 6 < 1 
mm. The width of the laser beam as it passes through the sample is 1-2 mm 
so changes in concentatlon over smaller distances than this are not 
detectable. Scans of absorption versus position in the liquid at various 
stages of growth were analyzed in an attempt to detect variations with time 
which would indicate slow convection, but no systematic charges could be 
seen and all our evidence suggests that mixing is rapid, presumably because 
of lateral temperature gradients. 


I 
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VI. DISTRI B UTION COEFFICIENT IN CZOCHRALSKI GROWTH 

Crystal growth by pulling from the melt allows the determination of k 
both by varying the pull rate V and the rotation rate m which 
determines the thickness of the boundary layer 6. 

Examples of crystals grown by pulling from the melt are shown In fig. 
15. The largest crystals grown were pulled from crucibles 4.2 cm In 
Internal diameter and In height and were up to 2.2 cm In diameter and lOOg 
In weight. Larger crystals could presumably be grown without difficulty by 
scaling up the apparatus. 

It is necessary to prevent thermal shock In order to avoid cracks and 
thus the pulled crystals were cooled very slowly and then annealed at 300’’C 
prior to fabrication. All the crystals have round shapes and no prominent 
facets. 

Three prominent directions were chosen to grow CsCdCla, which are 
[0001], [1010], and [10T2]. We have not had any major problems In growing 
It along any particular direction. For comparison, the [0001] and [1010] 
directions were chosen for this study of solute redistribution.. To observe 
the solute distribution change due to the velocity change effect only, the 
crystal diameter must be small enough not to be affected by any crucible 
wall effect; l.e., the melt is perfectly mixed by crystal rotation and 
solute convection in the melt near the crystal Interface Is not 
significantly influenced by the crucible wall. A rotation rate of 15 rpm 
was used for all crystals grown at different values of V, up to a maximum 
of 4 cm/s. The weight was between 5g and lOg from a total charge weight of 
50g. 

When crystals were grown at varying (o (fixed V), problems were 
encountered because of the change in diameter caused by changing convective 
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flow as the crystal rotation rate was varied. To study the rotation rate 
change effects on solute redistribution, we used rates from 2.2 rpm to 100 
rpm. At 120 rpm, the crystal exhibited an Irregular cross-section. 

The cobalt concentration C$ was measured by polishing the crystal 
along opposite sides and measuring the absorption of the He-Ne laser beam. 
The absorption was calibrated against samples analyzed by spark Sw ce 
emission spectroscopy which had been found to be accurate within ± 10%. 

The calibration gives the concentation in ppm as 220.4 a where a is the 
measured absorption coefficient. 

Data for <0001> and <10l0> orientations are shown in Figs. 16 and 17. 
The slope of Fig. 17 is 1.6 [3] so that can be 

evaluated using data for the kinematic viscosity v from Table 1. From 
the average slopes we obtain = 8.38 x 10“® cm^ s”^ from the <0001> 
data and Dj^ = 7.43 x 10“ ® cm^ s"^ for <1010> so the two values are in 
good agreement. 

To determine the values 8g/kT and Og for each orientation, we use 
the two sets of equations d/dV * -2396 sec/cm and 

An( (Yg+1 )/k^)^QQQ^^ = 1.92 for the <0001> orientation plus 
d/dV C^n(1/k-1)]^^Q;jQ^ = -1832 sec/cm and J^n(( Yg+1 )/k^ = 2.03 for 
the <10i0> orientation, all at V = 0.6 cm/hr. However, we do not have 
enough data to find the independent Dg value so we must use OgDg as a 
single value. Therefore, for the <0001> orientation 



Figure 16: Measured plots of fn(l/k-l> versus V for the Co - CsCdCl 
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From the above result, we know that the {0001} have the larger interface 
field strength which is what one might have expected because the {0001 } are 
the close packed planes in CsCdCls (H.C.P, structure). A summary of these 
determined properties is given in Table 3. 

It should be noted that distribution coefficient data for the same 

orientation are very different according to whether crystals are grown by 

gradient freezing or pulling from the melt. A direct comparison is shown 

in Fig. 18. The slopes are, by coincidence, very similar but what is 

striking is the large difference in intercept. In this case, with the data 

fitted in each case by a straight line chosen by the least squares method, 

the intercepts are k\j-Q = .027 for gradient freeze method and 

ky=0 = 0*136 for pulling from the melt. It is clear that the intercept 

cannot be the equilibrium value in both eases. An interpretation based on 

the interface field concept gives the intercept as ('^+Tg)/k^ with Yg an 

c 

incomplete gamma function and k^ “ k^[Yg(®)/Y5(0) 3 . Given that the 
distribution of solute ahead of the interface is different in the two 
cases, we might expect differences in the interface field to account for 
some of the discrepancy between the two intercepts. It is still doubtful, 
however, whether we can account for a factor of 5 when the system and 
orientation are the same. 

This large difference between the two sets of data is clearly an 
important result which shows that a single set of oistrlbutlon coefficient 
data cannot be relied upon to give the equilibrium value. Any detailed 
Interpretation is still speculative at the present time. 
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Table 3 

Summary of Physical Properties of CSCcIq ^^COq q^CI^ 
D|_ s (7.91 ± 0.48) X 10- ® cm^/sec 


> 

f 


k’= 0.115 
0 


< 0001 > 


< 1010 > 


Pg/kT -1.50 


-0.25 




3.98 X 10”^ cm/sec 


1.58 X 10"^ cm/sec 



coljalt in CsCd( 
CZ=Czo chralski 
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VII. SUMMARY AND CONCLUSIONS 

1 . The introduction of an interface field term due to local atomic 
configurations at the solid-liquid Interface has been shown to have a 
significant Influence on the theory of solute distribution. This 
effect can explain the observed curvature in the £n(k"^-1) versus 
growth rate data for germanium. 

2. The interface field effect has been introduced into the constitutional 
supercooling and morphological stability theories. Both stabilizing 
and destabilizing effects may result, and the magnitude of the 
resulting changes can be several orders of magnitude. 

3. A new model material for crystal growth studies, CsCdCls, has been 
grown by both gradient freeze and Czochralskl methods. Cobalt-doped 
melts in particular are very attractive for direct observational 
studies at temperatures where the infra-red radiation background is 
negligible. 

A. Studies of interface shape and the breakdown of stable growth have 
revealed a sensitivity of the interface curvature to the cobalt 
concentration in the liquid. This effect, attributed to radiation 
absorption in the liquid, has been, neglected in previous studies of 
this kind. 

5. A new technique for measurement of distribution coefficient by laser 
light absorption has been developed. 


6. Distribution coefficients during gradient freeze growth were found to 
be systamtically lower than those for Czochralski growth over the same 
range of growth rates. The values for the <0001> direction 
extrapolated to V = 0 differed by a factor 5, indicating a large 
difference from equilibrium in one or both cases. 
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